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1.  Objective 


The  objective  of  this  project  is  to  perform  analyses  and  designs  of  efficient  fault-tolerant 
quantum  circuits  for  incorporation  in  both  small-  and  large-scale  quantum  computers 


2.  Approach 


Although  large-scale  quantum  computers  do  not  yet  exist,  future  quantum  computers  have  great 
potential  for  solving  problems  beyond  the  capabilities  of  classical  computers.  However,  no 
general  methods  for  constructing  new  quantum  algorithms  are  known,  and  potential  quantum 
computer  showstoppers  presently  include  only  the  capabilities  of  factoring  huge  numbers  and 
finding  a  hidden  datum  in  a  huge  database.  Much  remains  to  be  discovered  concerning  the 
potential  power  of  quantum  computers. 

Any  quantum  computation  can  be  ideally  represented  by  a  unitary  transfonnation  acting  in  the 
Hilbert  space  of  the  computational  degrees  of  freedom  of  the  quantum  computer,  and  any  unitary 
transformation  can  be  faithfully  represented  by  a  network  of  universal  quantum  gates,  such  as 
two-qubit  controlled-NOT  (CNOT)  gates  and  single-qubit  gates.  This  is  the  basis  of  the 
quantum  circuit  model  of  quantum  computation  ( 1 ).  The  quantum  circuit  model  is  the  most 
widely  used  model  of  quantum  computation.  In  place  of  the  bits  of  classical  computers,  which 
assign  the  values  0  or  1,  the  quantum  circuit  model  is  based  on  qubits,  which  are  quantum 
systems  that  can  be  in  a  state  of  0  or  1  or  a  quantum  superposition  of  the  two,  and  can  be 
entangled  with  each  other.  Qubits  may  be  implemented  in  a  number  of  possible  ways:  photon 
polarization  states,  atomic  states,  quantum-dot  electron  states,  Josephson  junction  charge  or  flux 
states,  etc.  In  a  quantum  computer,  the  qubits  are  manipulated  by  a  network  of  quantum  gates  in 
such  a  way  that  their  final  state  corresponds  to  the  solution  to  a  computational  problem.  The 
quantum  gates  act  on  the  qubits  and  can  be  mathematically  represented  by  tensor  products  of 
Pauli  matrices. 

An  important  measure  of  quantum  circuit  complexity  and  the  difficulty  of  performing  a  quantum 
computation  is  the  number  of  quantum  gates  needed.  A  quantum  algorithm  is  considered 
efficient  if  the  number  of  required  gates  scales  only  polynomially  (not  exponentially)  with  the 
size  of  the  problem.  Quantum  circuit  networks  are  usually  analyzed  using  discrete  methods; 
however,  potentially  powerful  continuous  differential  geometric  methods  are  under  development, 
using  Riemannian  geometry.  Since  unitary  transformations  are  themselves  continuous,  this  is 
perhaps  not  a  surprising  development.  Using  these  differential  geometric  methods,  optimal  paths 
in  Hilbert  space  may  be  found  for  executing  a  quantum  computation  with  minimum  gate 
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operations.  An  appropriate  metric,  connection,  covariant  derivative,  curvature,  and  geodesic 
equation  must  be  formulated. 

A  new  innovative  differential  geometric  approach  to  quantum  circuit  complexity  was  recently 
introduced  by  Nielsen  et  al.  (2).  A  Riemannian  metric  was  fonnulated  on  the  space  of  multi¬ 
qubit  unitary  transformations,  such  that  the  metric  distance  between  the  identity  and  the  desired 
unitary  operator,  representing  the  quantum  computation,  is  equivalent  to  the  number  of  quantum 
gates  needed  to  represent  that  unitary  operator,  thereby  providing  a  measure  of  the  complexity 
associated  with  the  corresponding  quantum  computation.  The  Riemannian  metric  can  be  defined 
as  a  positive-definite  bilinear  form  defined  in  tenns  of  the  multi-qubit  Hamiltonian.  The  analytic 
form  of  the  metric  can  be  chosen  to  penalize  all  directions  on  the  manifold  not  easily  simulated 
by  local  gates.  In  this  way,  basic  differential  geometric  concepts  such  as  the  Levi-Civita 
connection,  geodesic  path,  Riemannian  curvature,  and  geodesic  equation  can  be  associated  with 
quantum  computation.  In  accord  with  the  Schrodinger  equation,  the  unitary  transformation 
expressing  the  quantum  evolution  is  an  exponential  involving  the  Hamiltonian.  The  Hamiltonian 
can  be  expressed  in  terms  of  tensor  products  of  the  Pauli  matrices,  which  act  on  the  qubits.  The 
geodesic  equation  in  the  manifold  follows  from  the  connection  and  detennines  the  local  optimal 
Hamiltonian  evolution  corresponding  to  the  unitary  transformation  representing  the  desired 
quantum  computation  and  minimizing  the  quantum  circuit  complexity. 


3.  Results 


A  Riemannian  metric  is  first  chosen  on  the  manifold  of  the  Lie  Group  SU(2n)  (special  unitary 
group  in  2n  dimensions)  of  //-qubit  unitary  operators  with  unit  determinant  (3-5).  A  traceless 
Hamiltonian  serves  as  a  tangent  vector  to  a  point  on  the  group  manifold  of  the  /7-qubit  unitary 
transformation  U.  The  Hamiltonian  H  is  an  element  of  the  Lie  algebra  ,s  w(2")  of  traceless  2”x2" 
Hermitian  matrices  (3),  and  is  taken  to  be  tangent  to  the  quantum  evolutionary  curve  c  lHtXJ  at 
t  =  0.  (Here  and  throughout,  I  chose  units  such  that  Planck’s  constant  divided  by  2n  is  ti  =  1 .) 

Independent  of  U,  the  Riemannian  metric  (inner  product),  ,  is  taken  to  be  a  positive  definite 
bilinear  fonn  ( H,J )  defined  on  tangent  vectors  (Hamiltonians)  H  and  J.  Following  (2),  the  /?- 

qubit  Hamiltonian  H  can  be  divided  into  two  parts  P{H)  and  Q(H),  where  P(H)  contains  only  one 
and  two-body  tenns,  and  Q(H)  contains  more  than  two-body  terms.  Thus 

H  =  P(H)  +  Q(H),  (1) 

in  which  P  and  Q  are  superoperators  (matrices)  acting  on  H,  and  obey  the  following  relations: 

P  +  Q  =  I,  PQ  =  QP  =  0 ,  P2=P,  Q2=Q.  (2) 

For  example,  in  the  case  of  a  three-qubit  Hamiltonian,  for  Pauli  matrices  c,  ,cr2,  and  cr3 ,  one  has 
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P(H)  =  x,  cr,  ®  /  ®  /  +  x2cr,  ®  /  ®  /  +  x3cr3  ®  /  ®  /  +  x4  /  ®  cr,  ®  / 

+  x5  /  ®  cr,  ®  /  +  x6  /  ®  cr3  ®  /  +  x7  /  ®  /  ®  cr,  +  x8  /  ®/  ®cr2 
+  x9  /  ®  /  ®  cr3  +  x10  <Tj  ®  cr2  ®  /  +  x,,  <Tj  ®  /  ®  cr2  +  x12  /  ®  Cj  ®  cr2 

+  x,3  cr2  ®  cr,  ®  /  +  x,4  cr2  ®  /  ®  <Tj  +  x15  /  ®  cr2  ®  <Tj  +  x,6  CTj  ®  cr3  ®  / 

+  x,7  cr,  ®  /  ®  cr3  +  x,8  /  ®  O'!  ®  cr3  +  x19  cr3  ®  CTj  ®  /  +  x20  cr3  ®  /  ®  cr, 

+  x21  /  ®  cr3  ®  (Tj  +  x22  <j2  <E)  a2  <E)  I  +  x22  a2  <E)  I  <E)  +  x24  /  ®  cr,  ®  cr3 

+  x25  cr3  ®  cr,  ®  /  +  x26  cr3  ®  /  ®  cr,  +  x27  /  ®  cr3  ®  cr,  +  x28  cr,  ®  cr,  ®  / 

+  x29  cr2  ®  cr,  ®  /  +  x30  cr,  ®  cr3  ®  /  +  x31  CTj  ®  /  ®  CTj  +  x32  cr2  ®  /  ®  cr2 

+  x33  cr3  ®  /  ®  cr3  +  x34  /  ®  Cj  ®  Cj  +  x35  I®ct2(8)<j2+Xti6  I  ®  cr3  ®  cr3 ,  (3) 

in  which  <8>  denotes  the  tensor  product,  and 


=  .r 
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Here,  all  possible  tensor  products  of  one  and  two-qubit  Pauli  matrix  operators  on  three  qubits 
appear  in  P{H),  and  analogously,  all  possible  tensor  products  of  three-qubit  operators  appear  in 
(Q(II).  I  exclude  the  tensor  products  including  only  the  identity  because  the  Hamiltonian  is  taken 
to  be  traceless.  Each  of  the  terms  in  equations  3  and  4  is  an  8x8  matrix.  The  various  tensor 
products  of  Pauli  matrices,  such  as  those  appearing  in  equations  3  and  4,  are  referred  to  as 
generalized  Pauli  matrices.  In  the  case  of  an  n-qubit  Hamiltonian,  there  are  4n  -  1  possible  tensor 
products  and  each  tenn  is  a  2nx2n  matrix. 


The  right-invariant  Riemannian  metric  for  tangent  vectors  H  and  J  is  given  by  (2) 

1 


H,J)  =  —Tr[HG(J)\,  G  =  P  +  qQ. 


(5) 
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Here  the  superoperator  G  is  defined  in  tenns  of  superoperators  P  and  Q,  and  <7  is  a  large  penalty 
parameter  which  taxes  more  than  two-body  terms.  The  length  /  of  an  evolutionary  path  on  the 
SU(2n)  manifold  is  given  by  the  integral  over  time  t  from  an  initial  time  tj  to  a  final  time  ty 
namely, 

l  =  \dt(H(t),H{t))V\  (6) 

and  is  a  measure  of  the  cost  of  applying  a  control  Hamiltonian  H(t )  along  the  path. 

Using  the  metric,  equation  5, 1  derived  from  first  principles  and  in  complete  detail  ( 6 ,  7)  the 
reduced  Riemannian  metric,  inverse  metric,  the  Levi-Civita  connection, 

K  =^Tr{G-l(p)([cr,G(T)]  +  [T,G(<Tm,  (?) 

and  the  following  expression  for  the  covariant  derivative  of  a  vector  Z  =  zCTcr  along  a  vector 

7(7  _  . 

=  y  a: 


VrZ  =  y°dZ/dx°+^{[Y,Z]  +  G1([Y,G(Z)]  +  [Z,G(Y)])},  (8) 

in  which  the  Einstein  sum  convention  is  understood.  In  equation  7,  Greek  letters  p  ,  o ,  and  r 
denote  generalized  Pauli  matrices  appearing  in  the  n-qubit  Hamiltonian,  and  the  notation  is  such 
that  a  Greek  index  (for  example,  p  in  T£t)  refers  to  a  particular  one  of  the  4"  - 1  possible 

generalized  Pauli  matrices.  An  important  step  in  my  derivation  was  to  obtain  the  following 
identity: 

X crTr(G  1  (<r)[r,  G(Z)])  =  2"  (T1  ([r,  G(A)]) .  (9) 


Next  considering  a  curve  passing  through  the  origin  with  tangent  vector  Y,  such  that  yc 


dxb 

dt 


it  follows  from  equation  8  that  the  covariant  derivative  along  the  curve  in  the  Hamiltonian 
representation  is  given  by 


DtZ  =VYZ  = 


-j-  +  yr{[Y,Z]  +  G1  ([Y,  G(Z)]  +  [Z,  G(Y)] } . 
dt  2 


(10) 


Because  of  the  right-invariance  of  the  metric,  equation  10  is  true  on  the  entire  manifold.  One  can 
next  proceed  to  obtain  the  geodesic  equation.  A  geodesic  in  SU (2" )  is  a  curve  U ( t )  with  tangent 
vector  H(t)  parallel  transported  along  the  curve,  namely, 

D,H  =  0.  (11) 
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However,  according  to  equation  10,  one  has 


DtH  =  +  i  {[H,H\  +  G  1  ([ H ,  G(H)}  +  [H,  G(H )])} , 

dt  2 

which,  when  substituting  equation  11,  becomes 


Next  defining 


L  =  G(H)  =  F  1  (H) , 


one  has 

^  =  2-G(ff)  =  of— 1 

dt  dt  \  dt  J 


Thus  substituting  equations  13  and  14  in  equation  15,  one  obtains 

d^  =  i[L,F(L)\. 
dt 


Equation  14  implies 


H  =  G\L), 


(12) 

(13) 

(14) 

(15) 

(16) 

(17) 


and  therefore  solving  equation  16  for  L  yields  the  Hamiltonian  H  producing  the  geodesic  path. 
Equation  16  is  the  geodesic  equation  for  the  locally  optimal  quantum  computational  paths  in  the 
SU (2”)  manifold.  It  is  a  nonlinear  first-order  differential  matrix  equation  of  the  same  form  as 
the  Lax  equation  for  Lax  pairs  L  and  iF(L).  Here  FI  and  L  are  2"  x  2"  matrices,  and  the 
superoperator  F  =  P  +  q  lQ . 

For  the  three-qubit  case  with  P  =  sS  +  T ,  in  which  S  and  T  are  superoperators  projecting  out 
one-  and  two-body  interactions,  respectively,  and  s  is  a  parameter,  the  solution  to  the  locally 
optimal  Hamiltonian  corresponding  to  the  solution  of  the  geodesic  equation,  equations  16  and 
17,  is  given  by 

H(t)  =  "1)soe''<1-0“I)OV£?<6j’  e-ini-P)(so+Qo)e-it(q  '-s  ')S‘>  +  cflei,("q '^SaQ  e  i,(q  ' 1)5°  (18) 

where  S0 ,  T0,  and  Q0  are  the  one-,  two-,  and  three-body  parts  of  the  initial  Hamiltonian. 


Also,  beginning  with  equation  7, 1  derived  in  complete  detail  (7,  8)  the  Riemannian  curvature 
and  the  sectional  curvature  of  the  SU(2")  group  manifold.  These  quantities  are  essential  for 
developing  increased  understanding  of  the  globally  optimum  quantum  computational  paths  in  the 
SU{2n)  group  manifold. 
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Other  work  performed,  but  currently  of  secondary  interest  and  which  because  of  page  limitations 
is  not  included  in  this  report,  includes  the  following: 

1 .  A  brief  survey  of  the  current  progress  in  developing  physical  implementations  of  quantum 
computers  and  identifying  the  issues  peculiar  to  each; 

2.  Mathematical  analyses  of  quantum  circuits  for  implementing  the  quantum  Fourier 
transfonn,  expressed  in  terms  of  matrix  algebra  associated  with  successive  stages  of  the 
algorithm; 

3.  Addressing  issues  of  fault  tolerance  and  the  associated  requirements  for  supplementary 
quantum  error  correction  circuits;  and 

4.  Detennination  and  analysis  of  possible  quantum  circuits  for  producing  generic  states  of 
quantum  entanglement. 


4.  Conclusions 


In  this  effort,  I  have  derived  the  necessary  differential  geometric  structure  of  the  577(2")  group 
manifold,  including  the  connection,  curvatures,  and  geodesic  equation  in  complete  detail  from 
first  principles.  This  fundamental  understanding  is  essential  for  the  proposed  practical 
applications.  It  is  a  necessary  step  in  the  design  of  optimal  quantum  circuits  for  implementing 
small-  and  large-scale  quantum  computation.  Solutions  to  the  geodesic  equation  need  extensive 
investigation,  keeping  in  mind  that  one  must  incorporate  piecewise-smooth  joining  of  local 
geodesics  in  order  to  produce  possible  global  geodesics  in  the  group  manifold.  Possible 
conjugate  points  must  be  incorporated  in  the  analysis,  and  possible  classical  computational 
obstructions  may  need  to  be  circumvented  in  order  to  facilitate  the  construction  of  globally 
optimum  geodesics  and  associated  minimum-complexity  circuits  for  large-scale  quantum 
computation. 
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6.  Transitions 


1 .  This  work  will  supplement  the  small  U.S.  Anny  Research  Laboratory  (ARL)-Sensors  and 
Electron  Devices  Directorate  (SEDD)  Mission  Program  on  Quantum  Computing. 

2.  I  have  initiated  collaborative  research  with  Dr.  John  Myers  and  Prof.  Tai  Tsun  Wu  of 
Harvard  University. 

3.  Five  papers  have  been  accepted  for  publication  (four  of  which  are  refereed). 

4.  I  presented  an  American  Mathematical  Society  (AMS)  Short  Course  lecture  on  this  work  at 
AMS  Meeting  in  Washington,  DC. 

5.  I  presented  work  at  the  Gordon  Research  Conference  on  quantum  computation. 

6.  I  gave  four  invited  talks  at  international  meetings. 

7.  I  periodically  informed  a  classified  Government  group  of  my  progress  on  this  research. 
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